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1.  Introduction 

The  evolution  of  discontinuities  in  solutions  of  nonlinear 
hyperbolic  equations  possessing  smooth  Initial  data  was  first  examined 
in  a  simple  problem  by  Riamann  [1].  His  conjecture'  concerned  the 
conditions  for  a  simple  wave  to  develop  a  discontinuity.  Ludford 
[2]  re-examined  this  conjecture  in  the  context  of  the  initial  value 
problem  for  unsteady  iaentropic  flow  of  a  perfect  gas.  This  unfoldinR 
process,  in  the  hodograph  plane,  provides  an  asymptotic  estimate  of 
the  time  to  breakdown  of  the  solution.  Zabusky  [3]  employed  the 
unfolding  method  to  determine  an  estimate  of  the  time  to  breakdown 
for  the  transverse  oecillations  of  a  nonlinear  model  string.  The 
work  of  Lax  [4]  and  Jeffrey  [5]  also  employs  the  Riemann  invariants  to 
develop  comparison  theorems  which  provide  upper  and  lower  bounds  for 
the  critical  time  of  singularity  occurrence.  An  alternate  method, 
characterized  by  its  simplicity  has  been  introduced  by  Ames  [6].  Large 
classes  of  quasilinear  equations  can  be  obtained  by  differentiation 
of  first  order  equations.  The  general  solutions  of  these  are  also 
solutions  of  the  corresponding  second  order  equations.  These  solutions 
display  a  critical  time  of  singularity  occurrence. 


This  critical  tins  analysis  is  applied  hare in  to  nonlinsar  wavs 


equations  which  result  from  rubber-like  materials  characterised  by 
Mooney-Rivlin  and  Neo-Hookean  bodies.  The  times  to  discontinuity 
evolution  from  smooth  initial  data  are  computed  and  used  to  ascertain 
the  region  of  validity  of  the  generalised  Lagrange  series  solution. 

2.  Fundamental  Equations 

The  present  investigation  pertains  to  straight  bars  or  wires 
with  negligible  transverse  dimensions  and  possessing  a  uniform  cross 
section  of  finite  area.  In  agreement  with  Nowinski  [7]  we  make  the 
following  additional  assumptions! 

(i)  Transverse  inertia  during  the  bar  motion  is  neglected. 

(ii)  In  compression  and  tension  zones  the  bar  does  not  experience 
material  instability, 

(iii)  The  material  is  perfectly  elastic  and  incompressible. 

(iv)  The  bar  is  subjected  to  simple  uni-directional  strain  in 
the  sense  that  the  only  identically  nonvanishing  stress 
component  is  the  longitudinal  normal  stress  component  which 
is  uniformly  distributed  over  the  cross  section, 

(v)  The  effect  of  strain-rate  on  the  constitutive  equations  is 
neglected  and  the  static  stress-strain  relations  are  extended 
to  the  dynamic  case. 

(vi)  The  bar  is  infinitely  long  so  that  no  reflections  of  waves 
occur  and  other  possible  wave  interferences  ars  discarded. 


Adopting  the  Lagrangian  formulation  lat  both  tha  matarial 
coord inata  X  and  spatial  coordinata  x  ba  rafarrad  to  tha  same  fixed 
Cartaaian  aystam,  one  of  whose  axes  coincides  with  the  axis  of  the 
bar.  Let  pQ  and  p  ba  tha  mass  densities  in  tha  stress  free  configur¬ 
ation  (associated  with  tha  X  coordinate)  and  deformed  configuration 
(associated  with  x  *  x  (Xtt)  ).  If  t  is  time,  oq  tha  normal  stress 
rafarrad  to  tha  undeformed  cross  section  of  tho  rod,  and  u  the  particle 
displacement,  than  Cauchy's  law  of  motion  becomes  (dames  and  Guth  [8]) 


bOp  _  j>v_ 
il  ~  it* 


Since 


V  _  hx 

lx 


X  =  I 


the  stretch  (extension  ratio) 
can  be  written  as 


Consequently,  Eq.  (1)  becomes 


Tt*  "  C 

or  in  terms  of  the  stretch, 


where 


_L  <l<% 

Po  d\  ) 

(we  assume 


cr,  =  cr0  ( x) 

ci<^/AA  >  O 
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Equations  (3)  and  (4)  ara  special  casaa  of  tha  ganaral  forme  traatad 
by  Ames  [6j. 

Undar  tha  assumption  that  tha  strain  energy  exists  the  theory  of 
finite  elastic  deformations  (Truasdall  [8],  Eq,  42.11)  furnishes  the 
strass-s tretch  relation 


for  an  incompressible  body  in  simple  extension.  Here  W  is  the 
elastic  strain  energy  function  and  for  an  incompressible  material  the 
strain  invariants  are 

i-2X"‘+xz;  zx  +  x‘2,  m '  i.  (v) 

From  the  experimental  results  of  Rivlin  and  Saunders  (see  Trueadell 
[8],  p . ( }14 )  experimental  data  is  wall  approximated  by 

W  =  «(l-3)  f(I  -3)  (8) 

where  a  is  a  constant  and  f  is  an  arbitrary  function  to  be  obtained. 

The  expanded  form 

W  =  V  1  {« 

with  6^  constants,  has  bean  employed.  Retention  of  only  the  linear 
term  leads  to  (10) 

W  =  0<U-3)  +  ^(E- l)  J  {«,{$><>) 


corresponding  to  the  Moonay-Rivlin  material.  If  B«6  we  obtain  Rivlin’s 
Neo -Hooke an  material. 


If  W  takas  tha  ganaral  form,  Eq.  (8),  than  Eq.  (4)  becomes 


at* 


teh'^x'4)  ♦ 


Mi:i*  ,!Mu  •sAellbc. 


Alternately,  we  may  write 


For  the  Kooney-Rivlin  material  Eq,  (12)  becomes 

And  for  the  Neo-Hookean  material,  Eq.  (12)  becomes 


(13) 


(14) 


In  Eq.  (14)  2o  has  been  replaced  by  E/3  a  value  which  is  suggested  by 
the  desireability  of  obtaining  the  familiar  infinitesimal  strain 
relation  =  from  Eq,  (6). 

From  his  own  theory  of  finite  elasticity  Seth  [9]  obtains  the 


corresponding  equation 


(15) 


which  differs  fundamentally  from  Eq.  (!»»). 

3,  Reduction  to  First  Order  Equations 

Equation  (11)  and  its  specializations  are  of  the  general  form 
given  by  Eq.  (4)  while  Eqs.  (12-15)  are  of  the  general  form  specified 
by  Eq.  (3).  As  previously  observed  both  are  special  cases  of  the  general 
quasilinear  equation 


(16) 
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treated  by  Ames  [6],  In  that  work  it  ia  shown  that  Eq.  (16)  raaulta  if 
ona  calculataa  r  and  s  derivatives  of  tha  general  firat  ordar  aquation 

K'S'*',  j--«*  <i?> 

and  eliminates  tha  cross  partial  darivativa  term  ur>.  Equations  of 
tha  fora 

<W  «*}„.=<>  (18) 

ara  obtainad  if  F=  »  and  -  $  (u)  .  In 

that  cast  tha  two  equations  for  F  become  F^,  ±4>(u)  fvj 
Clearly  F  s  P  are  solutions  for  F.  Upon  integrating  F=0, 

we  obtain  tha  general  solutions 

H  [uJ  ^(u;  -Cj 

and 

G-[u,  ft-  £(uj,dj  =0 

respectively,  where  H  and  G  ara  arbitrary. 

Alternatively,  aquations  of  tha  form 


(19) 

(20) 


uAA  ^ c0 


(21) 


are 


obtainad  if 

(the  assumption  (Fp/P^)*  -  generates  a  similar 

form  with  appearing).  Consequently,  tha  equations 

for  F  become  F>>  ±  <PW  z.Q  .  Solutions  for  F  are 
p  :  ^  p  S  <^/  and  upon  integrating  F*0  we  obtain 

h[un  ft+  A</>(u,j)  - o 


(22) 


and 


g[u„.;  A. ~  =•  O. 


(23) 


<*,  Calculation  of  Breakdown  "Tima" 


The  finite  "time"  to  tha  evolution  of  a  discontinuity  in  u^  or 
higher  order  derivatives  can  be  calculated  from  Eqs,  (19),  (20), 

(22)  and  (23),  For  example,  from  the  total  r  derivative  of  Eq,  (20) 
we  obtain 


“a 


G/i  +  _ 

Gw  ”  ajGu,  <&'(*)  } 


(24) 


which  becomes  arbitrarily  large  when 

Since  the  initial  "time"  of  this  occurrence  is  usually  of  most 
interest  we  write 


A*r  -  - %— 

Gu  $  (*■) 

where  the  minimum  is  estimated  over  the  appropriate  reage  of  the 
quantities  in  Eq.  (26),  (For  real  problems  we  are  interested  in 
positive  values  ).  On  occasion  the  general  solutions  are 

employable  in  simpler  forms.  For  example,  in  some  problems  the 
general  solution  may  be  used  as 
U  =  A  [a  + 

instead  of  Eq.  (19).  In  this  case  the  r  derivative  becomes 

A'(-) 


u*  = 


l- 


;  u)  c  A.  +  ^(w  )  As 


which  is  unbounded  when 


^  ~  /k\u>)  (p'iu) 


and  the  critical  (minimum) time  is 


"A  (u>J 


(25) 


(26) 


(27) 


(28) 


(29) 
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Similar  calculations  for  tha  aquation  U^^-0 

can  ba  carried  out  from  Eqs.  (22)  and  (23).  A  discontinuity  in 
the  first  derivative  can  be  discovered  by  inquiring  when  the 
second  derivative  becomes  unbounded.  From  Eq.  (23)  we  find 


(30) 


(31) 


implicit  nature  of  the  general  solutions  (Eqs.  (19),  (26),  (22)  and 
(23)  )  inhibits  their  use  in  the  determination  of  the  solutions. 
Alternately,  we  can  obtain  a  aeries  solution  to  the  equations  F>o 
which  are 

t  *>(u)  Ua -O  (32) 

in  the  first  case  discussed  and 

t  <^(Ua)  -O  03) 

in  the  second  case,  Lagrange  expansions  are  discussed  in  various 

contexts  in  Goursat  [10],  Bellman  [11],  Banta  [12]  (for  finite 
amplitude  sound  waves)  and  Ames  and  Jones  [13]  (for  a  Monge-Ampere 

equation). 

A  Lagrange  series  is  now  constructed  for  --0 

as  a  typical  example  of  the  methodology.  Suppose  u  has  a  Taylor's 
series  expansion  about  s=0, 


A> -o 


This  form  is  inconvenient  since  the  derivatives  are  with  respect 
to  s  and  not  r.  Replacement  of  the  a  derivatives  is  carried  out 
by  using  the  differential  equation  and  an  inductive  scheme  due 
to  Goursat  [10,  p,4G5]  (see  also  Banta  [12]).  If  UA  +$iu)  M*.®0 


then 


(35) 


Consequently,  from  Eqe.  (34)  and  (35),  for  /n.  % | 


(36) 


a  form  which  contains  only  r  derivatives  of  u.  If  the  process  of 
r  differentiation  and  evaluation  at  a«0  are  interchangeable  the 
series  takes  the  form 


where  U.(hyoJ  =  -f(*)  ,  the  "initial”  condition.  This 

series  is  valid  out  to  the  first  singularity  -  that  is  to  the 


smallest  breakdown  value  s  . 

c 


To  integrate  uo  +  ^Xur)  °  ©  from  ^q.  (33>  we  not®  that  ** 
becomes 

U*.)  =  O 

upon  differentiation  with  respect  to  r.  With  v  =  Eq.  (38) 


(37) 


(38) 


becomes 


rJ\,  t  —  O 


(39) 


an  aquation  taking  tha  same  form  previously  analysed.  Upon 
solving  for  v,  u  is  recovered  by  integrating  with  respect  to  r 


The  Lagrange  series,  Cq.  (37),  or  its  integrated  form 
describes  the  waveform  in  its  transition  from  the  smooth  initial 
form  f(r)  to  the  onset  of  breakdown, 

6.  Application  to  Rubberlike  Materials 

Results  of  the  preceding  analyses  are  applicable  to  a  wide 
variety  of  problems.  Herein  we  will  investigate  Cq.  (4)  as  it 
applies  to  Mooney-Rivl in  materials.  For  that  application  Cq.  (4) 
becomes 

From  Section  3  the  first  order  equations  generating  Cq.  (40)  aim 
obtained  by  setting  U.  -  X  ;  A  (  \sj 
and  <P  -  [-^r  [*(/«■  + -*/3  j  ^ 

whereupon 

X*  t  <p(\) -o, 

If  the  initial  stretch  is  provided  by  X(XiO)  =  f(£)  then 
the  general  solution  of  Eqs.  (41)  become  (see  Eq.  (27)  ) 

X  =  f[x  *£(X)t] 


(40) 


(41) 


(42) 


»o  ^  ’I 
calculated  from  the  positive  minimum  of 


the  breakdown  time  can  be 


if  o  <  Xrt  <  X  <  X 

n  the  posit i 

t  =  t (i ftp  ,-nWt 


(43) 


If  {/(id)  $  A  and  tha  material  is  Neo-Hookean  (  £  =  o) 

then  (recall  e(  =  E/u  ) 


if  \0  a  |  -i..e.the  material  is  always  in  tension. 

The  Lagrange  series  becomes,  in  this  case, 

x(v) 

which  can  be  employed  to  study  the  onset  of  the  discontinuity. 
Since  X  —  +  -|ijr  the  x  and  u 

are  recovered  by  integrations. 


References 


1.  Riemann,  B.,  "Uber  die  Fortpflanzung  Ebener  Luftwellen  von 
Endlicher  Schwingungswerte",  Abh.  Gas.  Wise.  Gdttingen  8;, 

43,  I860, 

2.  Ludford,  G.S.S.,  "On  an  Extension  of  Riamann'a  Method  of 
Integration  with  Applications  to  One-Dimensional  Gas  Dynamics", 
Proc.  Carab.  Phil.  Soc.  48.,  499,  19S2. 

3.  Zabusky,  N.  J,,  "The  Exact  Solution  for  the  Vibration  of  Nonlinear 
Continuous  Modal  String,"  J.  Math,  Phys,  3,  1028,  1962, 

4.  Lax,  P.  D.,  "Development  of  Singularities  of  Solutions  of  Nonlinear 
Hyperbolic  Partial  Differential  Equations,"  J,  Math.  Phys.  5^ 

611,  1964. 

5.  Jeffrey,  A.,  "The  Evolution  of  Discontinuities  in  Solutions  of 
Homogeneous  Nonlinear  Hyperbolic  Equations  having  Smooth  Initial 
Data,"  J.  Math.  Mech.  17,-331,  1967. 

6.  Ames,  W.  F.,  "Discontinuity  Formation  in  Solutions  of  Homogeneous 
Nonlinear  Hyperbolic  Equations  Possessing  Smooth  Initial  Dat a'i 
Submitted  to  U.S,  National  Congress  of  Applied  Mechanics. 

7.  Now inski,  J.  L*,  "On  the  Propagation  of  Finite  Disturbances 
in  Bars  of  Rubber-like  Materials",  ASME  Paper  64-WA/RP-9. 

8.  Truesdell,  C.,  "The  Mechanical  Foundations  of  Elasticity  and 
Fluid  Mechanics",  Jour.  Rat,  Mech,  Anal.  1,  125,  19S2. 


Seth,  B.  R, ,  "Finite  Longitudinal  Vibrations"*  Proc,  Indian 
Acad.  Sci.  A2S.  151,  1948. 

Goursat ,  E«,  "A  Coursa  in  Mathematical  Analysit",  Vol.  I, 
Dovar,  New  York,  1959. 

Bellman,  R.  E.,  "Perturbation  Technique*  in  Hathamatics, 
Physics  and  Engineering"  (Athena  Series),  Holt,  Rinehart 
and  Winston,  New  York,  1964. 

Banta,  E.  0.,  "Lossless  Propagation  of  One-Dimensional  Finite 
Amplitude  Sound  Waves",  J.  Math.  Anal.  Appl,  10,  166,  1965. 
Ames,  W,  F.  and  Jonas,  S,  E.,  "Integrated  Lagrange  Expansions 
for  a  Mongo-Ampere  Equation",  J»  Math.  Anal.  Appl.  21, 

479,  1968. 


Security  CU»«Iflc»tlon 


DOCUMENT  CONTROL  DATA  -RAD 

(Security  dtiifflcetew  of  ilflm.  body  of  ebanect  and  Inditing  annotation  muat  ba  an  tar  ad  atban  trie  qyreli  met  It  c 


I.  ORIGINATING  ACTIVITY  fCoipoMt,  author) 

The  UniverBity  of  Iowa 


U.  RKRORT  (KCURlTY  CLAMIRIC  ATION 

unclasBlfled 


Sb.  GROUP 

NA 


».  RIPORT  TITL* 


On  Wave  Propagation  in  One  Dimensional  Rubberlike  Materials 


o**CRiRTivr.  NOT**  f  7>pe  of  report  end  incfusfv*  dafmaj 

Technical  Report 


*■  autnoR(I)  ml  dot*  initial,  name) 

William  F,  Ames 


•  ■  RIPORT  OAT* 

October  1969 


7a.  TOTAL.  NO-  OF  PACK! 

13 


76.  NO.  OF  R*fr* 

lU 


M.  CONTRACT  Oft  GRANT  NO- 

DAAF03-69-C-001U 

b.  PROJECT  NO. 


»4.  ORIGINATOR**  RIPORT  NUhf»*R(*> 

Project  Thenis  Report  No.  lL 


*6.  othcr  riport  nO(I)  (Any  other  number*  Hmt  mmy  be  eaaigned 
thia  report) 


10  DISTRIBUTION  ITATEMPST 

Distribution  of  This  Report  is  Unlimite 

d. 

tt-  supplement  \  rv  note* 

13.  XRONAORING  MILITARY  ACTIVITY 

U.S.  Army  Weapons  Consnand 

RLE  Directorate 

Rock  Island,  Illinois  61201 

I..  abstract 


Two  methods  for  studying  the  evolution  of  dicontinuities  in 
solutions  of  nonlinear  hyperbolic  equations  possessing  smooth  initial 
data  are  extant.  A  Third  method,  characterized  by  its  simplicity, 
results  from  the  observation  that  large  classes  of  quasilinear 
equations  can  be  obtained  by  differentiation  of  first  order  nonlinear 
equations.  The  first  order  equations  are  employed  to  calculate  the 
time  to  singularity  occurrence.  This  result  is  generalized  and 
applied  herein  to  nonlinear  wave  equations  from  one-dimensional  rubber¬ 
like  materials.  The  time  to  discontinuity  evolution  from  smooth 
initial  data  is  computed  and  used  to  ascertain  the  region  of  validity 
for  a  generalized  Lagrange  series  solution. 
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